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ABSTRACT 

We compute the expected luminosity function of GRBs in the context of the inter- 
nal shock model. We assume that GRB central engines generate relativistic outflows 
characterized by the respective distributions of injected kinetic power E and contrast 
in Lorentz factor k — rmax/rmin- We find that if the distribution of contrast extends 
down to values close to unity (i.e. if both highly variable and smooth outflows can 
exist) the luminosity function has two branches. At high luminosity it follows the dis- 
tribution of E while at low luminosity it is close to a power law of slope —0.5. We then 
examine if existing data can constrain the luminosity function. Using the log N — log P 
curve, the E^ distribution of bright B ATSE bursts and the XRF / GRB ratio obtained 
by HETE2 we show that single and broken power-laws can provide equally good fits 
of these data. Present observations are therefore unable to favor one form of the other. 
However when a broken power-law is adopted they clearly indicate a low luminosity 
slope ~ —0.6 ± 0.2, compatible with the prediction of the internal shock model. 
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1 INTRODUCTION 

The isotropic luminosity of long gamma-ray bursts is known 
to cover a wide range from underluminous, nearby bursts 
such as GRB 980425 or GRB 060218 (with L < lO"*^ 
erg.s"^) to ultrabright objects like GRB 990123 (L > 10^^ 
erg.s"^). While it has been suggested that the weak est GRBs 
could simply be normal events seen off-axis (jYamazaki et al.l 
l2003l ). this possibility has been recently discarded both from 
limits on afterglow br i ghtness and for statistical r easons 
l|Soderberg et al.ll2004l : iPaigne fc MochkovitchI |2007| ). The 
difference of six orders of magnitude between the brightest 
and weakest GRBs is therefore probably real. The parame- 
ters (stellar rotation, metallicity, etc.) which are responsible 
for this diversity in radiated power are not known. However, 
in the restricted range lO'^^ <L< 10^^ erg.s"^ the value of 
the isotropic luminosity is possibly fixed by the opening an- 
gle of t he jet which may a lways carry the same characteristic 
energy dMTeEaLlllool). 

The purpose of this paper is to see how basic theoretical 
ideas and existing data can be used to constrain the GRB lu- 
minosity function (hereafter LF) p{L). First, we should insist 
that p{L) here represents the "apparent" LF which includes 
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viewing angle effects and beaming statistics (i.e. bursts with 
narrow jets are more likely seen off-axis and therefore under- 
represented in the distribution) . It is therefore different from 
the "intrinsic" LF po{L) which would be obtained with all 
GRBs seen on-axis. In the lack of a complete, volume limited 
sample of GRBs with known redshift, only indirect observa- 
tional indicators such as the log TV — log P plot can constrain 
the LF. These indicators however depend not only on p{L) 
but also on the GRB rate and spectral energy distribution. 
The simplest possible form for p{L) is a single power law 
p{L) oc between Lmin and I/max. Together with the pa- 
rameters describing the GRB rate and spectral shape, S, 
Lmin and Z/max cau be adjusted to provide the best possible 
fit of the available indicators. Considering the mixing of the 
LF with other quantities in the fitting process it is remark- 
able that studies using different observational constraints 
have converged to a similar value of the slope 5 ^ 1.5 - 1.7 
(e.g. iFirmani et "al]l2004l : iDaigne fc Mochkovitchll2007l ). 

In a second step one can consider the more general case 
of a broken power law LF with now five parameters: Lmin, 
Lmax, the two slopes 5i and 82 and the break luminosity 
Lb. We will see in Sect. 2 that there is some indication that 
the internal shock model of GRBs can produce a broken 
power law LF and we want to check if it is also favored by 
the existing observational data. As in our previous study we 
have used a Monte Carlo method to generate a large number 
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of synthetic events where the parameters defining the burst 
properties are varied within fixed intervals. Preferred values 
of the parameters are those which yield the minimum 
for a given set of observational constraints. We summarize 
these constraints and present the Monte Carlo simulations 
in Sect. 3. We discuss our results in Sect. 4 and Sect. 5 is our 
conclusion. 



For Lmin ^ -£/ ^ Lmax, the probability to have a radiated 
luminosity in the interval [L; L + dL] is po{L)dL, where the 
intrinsic LF po{L) is given by 



Po (L) = 



dK . (5) 



THE GRB LUMINOSITY FUNCTION IN 
THE INTERNAL SHOCK MODEL 



The internal shock model (|Rees fc Meszarod [l993) is the 
most discussed solution to explain the prompt gamma-ray 
emission in GRBs. In this section, we demonstrate that 
it naturally leads to a broken power-law LF with a low- 
luminosity slope close to —0.5. 



2.1 



Internal shock efficiency and luminosity 
function 



In th e context of the internal shock model (|Rees fc Meszaroi 
1 1994 ). the prompt gamma-ray emission is produced by rel- 
ativistic electrons accelerated in internal shocks propagat- 
ing within a relativistic outfiow. The (isotropic equivalent) 
radiated luminosity L is then a fraction of the (isotropic 
equivalent) rate of kinetic energy injected in the flow E 

L^f{^)E. (1) 

The efficiency / (k) is the product of three terms 

/ («:) - /rad /dyn (t) , (2) 

where (i) /dyn («:) is the fraction of the kinetic energy 
which is converted by internal shocks into internal energy 
("dynamical efficiency"). This fraction depends mainly 
on the contrast k = rmax/rmin of the Lorentz factor 
distribution in the relativistic outflow; (ii) to is the fraction 
of this dissipated energy which is injected into relativistic 
electrons. We assume that Ec is close to the equipartition 
value Ec = 1/3, as it is a necessary condition to have 
an efficient mechanism; {iit) /rad is the fraction of the 
electron energy which is radiated. To explain the observed 
variability timescales in GRB lightcurves and to maintain 
a reasonable efficiency, the relativistic electrons must be in 
the fast cooling regime, which means that their radiative 
timescale is very short compared to the hydrodynamical 
timescales in the outflow. In this case, we have /rad — 1- 

The GRB LF in the internal shock model is therefore 
related to the physics of the relativistic ejection by the cen- 
tral engine. We assume that the contrast n is distributed 
between ACmin and Kmax with a density of probability t/j (k) 
and that E is distributed between -Bmin and Emax with a 
density of probability (f) [E^ . The minimum and maximum 
radiated luminosities are therefore 



Lmin — / (/^min) -E'it 



and 



/(Kh 



(3) 



(4) 



2.2 The case of a power-law distribution of 
kinetic energy flux 

We assume that the injection rate of kinetic energy in the 
relativistic outflow follows a power-law distribution 



E) 



1 



\E mi 11 



E 



(6) 



with -Emax ^ ^'min. Thcu, the GRB LF given by equa- 
tion ((5]) becomes 

5-l(L\-' 
Po(L)-^(^j X 

where the luminosity L* is deflned by 

L, = / (Kmax) -Emin • (8) 

Let us now consider a flrst case where GRB central en- 
gines can produce all kinds of outflows, from highly vari- 
able to perfectly smooth. In this case, the minimum con- 
trast is Kmin = 1, corresponding to a minimum efficiency 
/ (Kmin) = 0, as uo internal shocks can develop in an outflow 
with a constant Lorentz factor. The first limit in the inte- 
gral in equation ((7)l is then always given by / (L/ijmax) 
and the LF is made of two branches : 

• High-luminosity branch. For L, ^ L ^ imax, the 
second limit in the integral is Kmax, which leads to 

Po [L] ~ (^^^ f (Kmax) X 



dx X 



'4>[ f {Xf (ftmax)) 



/' / (a;/(«max)) 



(9) 



For L* ^ L ^ imax, i/imax — >■ SO that the integral 
is nearly constant. The high-luminosity branch of the LF is 
therefore very close to a power-law of slope —S, i.e. the slope 
of the intrinsic distribution of injected kinetic power. 

• Low-luminosity branch. For imin = 0^1/^1/*, 

the second limit in the integral is / (i/iSmin), which leads 
to 



PO (L) 



J I'^maxJ 



L/L„ 



dx X 



V' / (a;/(Kmax)) 



(10) 
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The LF is determined by the behaviour of the efficiency / (k) 
at very low contrast. We assume that / (k) ~ /o (k — 1)" for 
K — >■ 1 and we define Lo = /o-Emin. Then, for L ^ L,, we 
have 



PO {L) 



1 



V.(i) 



l + a((5 - 1) 



(11) 



We therefore find that for L <^ and ^min = 1, the LF 
is also a power-law, with however a slope 1/a — 1 which 
does not depend on the slope of the intrinsic distribution of 
injected kinetic power. 

Then, the predicted intrinsic GRB LF in the internal 
shock model is a broken power-law, of slope 1/q — 1 at 
low-luminosity and —S at high luminosity, with a break 
luminosity L, = / (^max) i^min. The shape of po{L) at the 
transition {L ~ L,) is related to the distribution of the 
contrast ^(k). If fCmin 7^ 1 this result remains valid as long 
as very low contrasts can be achieved (fCmin ^ 1.1; see 
figure [U left panel). Note that the a nalysis of the interna l 
shock model parameter space made bv lBarraud et al.l (|2005l ') 
shows that very low contrasts are necessary to produce soft 
GRBs such as X-ray rich GRBs (XRRs) and X-ray hashes 
(XRFs). However, if it happens that GRB central engines 
never produce smooth outflows (i.e. if Kmin is not close to 
unity), the calculation made above remains valid for the 
high-luminosity branch, which is still a power-law of slope 
—5, but the low-luminosity branch is much reduced and no 
more a power-law. 

Finally, we have briefly considered the case where the 
distribution of injected kinetic power (j>{E) is not a power- 
law. For example, for a log-normal distribution peaking at 
Et,, we again obtain that the instrinsic LF follows <^(-E) 
at high luminosity and is a power-law of slope —0.5 at 
low-luminosity, with a transition at L, ~ / (fimax) £*• 



2.3 A simple model for the internal shock 
efficiency 

To investigate more precisely the GRB LF, we need 
to know the form of the e fficiency /(k). As shown in 
iDaigne Sz MochkovitchI l|l998h , it is a priori a complex func- 
tion of the initial distribution of the Lorentz factor and the 
kinetic energy in the relativistic outfiow. However, one can 
make a simple estimate by usi n g the toy mode l devel oped in 
iDaigne fc MochkovitchI (|2003l '): lBarraud et al.l (|2005l ) where 
we only consider direct collisions between two equal-mass 
relativistic shells. In this case, the efficiency is simply given 
by 



/(^) 



K+l 



(12) 



For low contrast, it behaves as / (k) ~ Ec (k — 1) /8, which 
corresponds to fo — ec/8 and q = 2. In addition to 
this toy model which gives an explicit expression for the 
efficiency we have used our detailed internal shock code 
l|Daigne Sz MochkovitchI Il998l) and checked the quadratic 
dependence of /{k) in {ft — 1). Also notice that the result 
a = 2 will remain valid even if is not strictly constant, as 
long as it does not vary as some power of (k — 1) at low k 



(te oc (k— 1) leading for example to a = 3 andpo(i) oc L"^''^ 
at low luminosity). 

Assuming a constant e^, we have plotted the result- 
ing GRB LF in figure [1] The intrinsic distribution of in- 
jected kinetic power is defined between -Bmin = 10^^ erg s~^ 
and -Emax ~ 10^* erg s~^ and has a slope —S = —1.7. We 
have fixed the maximum value of the contrast to ftmax = 
10, so that L« ~ 1.4 x 10^^ erg s~^. In the left panel, 
we have assumed that the logarithm of the contrast k is 
uiuformly distributed between log Kmin and log Kmax , with 
Kmin = 1, 1.001, 1.01, 1.1 Or 2. For all values of Kniin, the high- 
luminosity branch is the same power-law of slope —5 = —1.7. 
For Kniin ~ 1, the low- luminosity branch extends down to 
L = and is the expected power-law of slope 1/a — 1 = 
— 1/2. This branch is still clearly visible for Kmin ~ 1.001, 
1.01 and 1.1 but nearly disappears for Kmin = 2. For even 
higher values of Kmin, only the high- luminosity power- law 
remains. 

We have tested in the right panel of figure [T] that for other 
choices of the distribution of contrast ii{K), the GRB LF is 
not affected (the two slopes remain unchanged) except for 
the shape of the transition at E ~ L». A low-luminosity 
branch of slope ~ —0.5 in the intrinsic LF is therefore a ro- 
bust prediction of the internal shock model, as long as GRB 
central engines can produce smooth outflows (very low con- 
trasts). The low- luminosity branch will however manifest 
itself only if E, = /(Kmax)£min is large enough; otherwise 
the observationally accessible part of the LF of cosmolog- 
ical GRBs will behave as a single power-law (in figure [1] 
L* ~ 1.4 X lO'^^ erg s"^). 



2.4 Apparent GRB luminosity function 

The GRB LF that has been derived from the internal shock 
model is intrinsic. If GRB ejecta have a jet-like geometry 
with an opening angle AO which is not correlated to the 
kinetic energy fiux E, the apparent LF above E* has the 
same shape as the intrinsic one since the fraction of observed 
GRBs does not depend on E. At lower luminosities, two 
effects are in competition : low-luminosity bursts can be due 
to a low internal shock efficiency and/or a large viewing 
angle. Close to E«, the first effect dominate and the slope 
is still very close to —0.5 as predicted above. At very low 
luminosity, the second effect takes over. It can be shown that 
the slope then becomes close to —7/6. Observing this final 
slope seems difficult as it involves the detection of very faint 
bursts. However viewing angle effects already modify the 
low-luminosity slope below the break where it progressively 
departs from its intrinsic value —0.5 (see Fig. 2). 

If it happens that the opening angle A9 is correlated 
with E, the apparent p{L) will be also different from the 
intrinsic one at high luminosity. A possible correlation could 
be 



E(l - cos Ae^) = est = En 



(13) 



meaning that the true kinetic energy rate is the same for 
all bursts. Such an assumption is motivated by the evidence 
that there might be a standard energy reservoir in GRBs 
(|Frail et al.ll200l1 ). In this case, the high-luminosity branch 
(above E*) of p{L) has a slope — (1 -I- S), where —5 is the 
slope of the intrinsic LF. This is illustrated in figure [2] (right 
panel) . 
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Figure 1. The intrinsic GRB LF in the internal shock model. The function po{L) is plotted, assuming an intrinsic distribution 
of injected kinetic power E that is a power-law od slope —5 = —1.7 between i?niin = 10^^ erg and -Bmax = 10^* erg s~^. This 
function <^(-E) is plotted as a thin line in the left panel. Left: effect of the minimum constrast Kniin- We assume that log/t is uniformly 
distributed between logKmin ^nd and logftmax with ftmin = li 1.001, 1.01, 1.1 and 2, and Kmax = 10. Right panel: effect of the shape of 
the distribution of constrast. We fix Kmin = 1 and Kmax = 10 and consider three different shapes for the distribution of constrast ■>P{k) : 
(i) logK uniformly distributed (solid line); (ii) ^(k) oc exp(— k/4) (dotted line); (iii) ft uniformly distributed (dashed line). An inset 
shows the transition at L* = 1.4 X 10^^ erg in more details. 




Figure 2. The apparent GRB LF in the internal shock model. The intrinsic LF po{L) is plotted in dashed line for log(K) uniformly 
distributed with fCmin = 1 B.nd Kmax = 10, and for E following a power-law distribution of slope — <5 between Emin = 10^^ erg s~^ and 
i-max = 10^* erg s"'^. The corresponding apparent LF is plotted in solid line. Left: the opening angle is distributed between and n/2 
according to p{A6) = sin (uniformly distributed opening angle). We adopt S = —1.7. Right: the opening angle is correlated with the 
kinetic energy flux according to Eq. II13I I. We adopt S = —0.7. 
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3 CONSTRAINING THE GRB LUMINOSITY 
FUNCTION 

3.1 Monte Carlo simulations 

We use Monte Carlo simulations to constrain the GRB LF 
from observa t ions. The method is described in details in 
IDaigne et al] (|2006l ). We recall here the main lines : 

• Properties of the long GRB population. We char- 
acterize this population by the intrinsic distribution of sev- 
eral physical properties : (1) the comoving rate. We de- 
fine 7?.grb(2) (Mpc~^ yr~^) as the GRB comoving rate 
at redshift z. We assume that the GRB comoving rate is 
proportional to th e star formation rate (SFR). Following 
iPorciani &: Madaul (|200ll '). this can be written as 7?.grb = 
k X TZsN, where 7?.sn is the comoving rate of collapses of mas- 
sive stars a bove 8 . We cons ider three possible scenarios 
(see Fig.l in lDaigne et al.l l|2006f ) and reference therein), that 
all fit the observed SFR up to « ~ 2 — 3 : SFRi where the 
SFR decreases for z >2, SFR2 where it is constant for z >2 
and SFR3 where it increases for 2: > 2. In this last case, we 
have to assume a maximum redshift for star formation. We 
adopt Zniax = 20; (2) the LF. In this paper, for simplicity, we 
do not discuss evolutionary effects. Therefore, the probabil- 
ity density p{L) of t he isotropic equivale nt luminosity L does 
not depend on z. In lDaigne et al.l l|2006h we only considered 
the case of a power-law distribution, with p{L) oc for 
imin i ^ imax . Here, wc test more complicated LFs, i.e. 
broken power-laws defined by 



p(L) 



for Lniin ^ 1/ ^ -Lb 
for Lb 5? -L < I/max 



(14) 



(3) the distributions of intrinsic spectral parameters. We as- 
sume that the GRB photon spectru m is given by a br oken 
power-law with a break at energy Ep (|Band et al.lll993 ) and 
a slope —a (resp. — /?) at low (resp. high) energy. We checked 
th at the use of the m ore realistic spectrum shape proposed 
bv lBand et al] l| 19931 ) does not affect our conclusions. The 
slopes Q and P are given the distribu tions observed in a 
sample of BATSE bright long GRBs (|Preece et al.ll200(j ). 
For the peak energy Ep we consider two possible cases, ei- 
ther a log-normal distribution, with a mean value Ep^ and 
a dispersion a = 0.3 dex (hereafter "log-normal Ep distribu- 
tion") or an intrinsic correlation between the spectral prop- 
erties and t he luminosity (here after "Amati-like relation"), 
as found bv lAmati et all l|2002l '): lAmatil (|2006l 'l. We assume 
in this case that 



Ep = 380 keV 



1.6 X 1052 erg s-i 



(15) 



with a nor mal dispersion a = 0.2 d ex in agreement with ob - 
servations (|Yonetoku et al' ('2004'') : ' Chirlanda et al? ('20051), 



Band fc Pr cccc (2005) 



see however [Nakar & Piran (2005); 
who tested this relation against BATSE data and concluded 
that selection effects were dominant). 

• Criteria of detection by several instruments. 
With the assumptions listed above, a GRB in the simu- 
lation is characterized by a redshift 2, a peak luminosity 
L, and a spectrum defined by Ep, a and /?. It is therefore 
possible to compute the observed peak flux in any spectral 
band. Using the kno wn sensitivity of several instruments (see 
IDaigne et al.l l|2006l ) for the detailed thresholds we use), we 



can determine if a given burst is detected by the following 
experiments : {1} BATSE, for which we define two samples in 
our synthetic bursts (all BATSE bursts and bright BATSE 
bursts) ; (tt) HETE2, for which we test the detection ei- 
ther by the gamma-ray (FREGATE) or the X-ray (WXM) 
instrument; and (iii) SWIFT, for which we test the detec- 
tion by the gamma-ray instrument only (BAT) and we also 
define two samples (all SWIFT bursts and bright SWIFT 
bursts). It is then possible to compute simulated observed 
distribution of various quantities to compare them with real 
data. 

• Observational constrains. We use three different 
kind of observations: (1) the log A^— log P diagram of BATSE 
bursts l|Stern et aLll2000l . lioOll '): (2) th e observed peak en - 
ergy distribution of long bright GRBs l|Preece et al. liooO); 
and (3) the observed fraction of soft GRBs (X-ray Flashes 
and X-r ay rich GRBs) in the sample of GRBs detected by 
HETE2 (; Sakamoto et all 120051 ). The logA^ - logP diagram 
is broadly used for s uch kind of studies but we have shown in 
IDaigne et al.l (|2006l ) that (2) and (3) are good complements 
to better constrain the parameters of the GRB population. 

Depending on the assumptions on the spectral properties, 
we have 4 or 5 free parameters for a single power-law LF 
(fc, Lmin, imax, S and Epfi iu the case of a log-normal 
distribution). The observational constraints correspond 
to 41 data points (see figure U) : 30 data points in the 
log iV - log P diagram published by IStern et al"] (120011): 10 
points for the Ep distribution published bv IPreece et al] 
(|2000l ) which has been rebinned in 10 logarithmic bins of 
size 0.2 dex between 15.8 keV and 1.58 MeV; 1 point for the 
fraction of soft GRBs). We have therefore 37 or 36 degrees 
of freedom. The numerical procedure is the following : for 
a set of parameters, we generate randomly a population of 
10^ GRBs using the distribution defined above, we then 
compute the simulated distributions of observed peak fiux, 
peak energy, etc. and we compare them to real data, by 
computing a x^- We do that for a large number of sets of 
parameters, randomly chosen to explore a large space. We 
always find a clear minimum Xmin of a-nd we define as 
"best models" all models with Xmin ^ ^ Xmin + ^X^j 
where Ax^ defines the la level and is computed from the 
number of degrees of freedom. 



The main results obtained in IDaigne et al.l l|2006l ) are 

(1) that SFR3 is strongly favored by the observed redshift 
distribution of SWIFT bursts. But a SFR rising at large z 
appears unlikely as it would overproduce metals at early 
times. This is therefore an indirect indication in favor of 
a GRB rate that does not directly follow the SFR, for 
instance due to an evolution with redshift of the efficiency 
of massive stars to produce GRB^; (2) with this SFR3, 
both the "Amati-like relation" and the "log-normal Ep 
distribution" give good fits to the observational constraints 



^ In jPaigne et al. I 1 I2OO6I) . we have tested whether an evolution 
of the LF could reconcile the SFRi or SFR2 scenario with Swift 
data. We found that this is very unlikely as the evolution should 
be strong (L oc (1 -|- z)" with u > 2). We therefore conclude that 
the evolution of the GRB rate (i.e the evolution of the stellar 
efficiency to produce GRBs) is dominant compared to a possible 
evolution of the LF. 
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listed above. Best model parameters of the LF do not 
vary too much from one scenario to another : the slope 
is well determined, 5 ~ 1.5 — 1.7, but the minimum and 
maximum luminosities are not so well constrained, with 
Lmin ~ 0.8-3x10™ erg s"^ and L^ax ~ 3-5x10^^^ erg s"^ 
It is interesting to note that with a different methodology, 
several groups have confirmed our conclusions on the GRB 
comoving ra.te jLe fc DermeJ l2007l : iGuetta fc PiranI l2007l : 
iKistler et al.ll2007l ). 

In this paper, we present the results of additional sim- 
ulations that we have carried out to test if the GRB LF can 
be a broken power-law. With two supplementary parameters 
(the break luminosity Lb and the second slope), we have now 
6 (resp. 7) free parameters in the Amati-like relation case 
(resp. the case of a log-normal peak energy distribution). It 
is difficult to constrain accurately so many parameters with 
Monte Carlo simulations. Therefore, we have chosen to keep 
the maximum luminosity constant in all our simulations. 
We ad opt Lmax = 10^^'^ erg s~^, according to our previous 
study ijDaigne et al.|[200^ ). We also keep Lmin constant, and 
equal to a low value corresponding to weak GRBs that can- 
not be detected at cosmological distance. We usually adopt 
10*^ erg s-i but we have also tested other values 
(see next section). Keeping these two luminosities constant 
in our Monte Carlo simulation, we have the same number of 
degrees of freedom than in the model with a simple power- 
law LF. 

3.2 Results 

In our whole new set of simulations, we always find a clear 
minimum of . The parameters of the best model, as well 
as 1 cr error bars are listed in table [T] As can be seen, we 
focus on the scenario where the comoving GRB rate follows 
SFR3 and the peak energy is given by the Amati-like re- 
lation. For comparison, we also give the para meters of two 
refere nce models with a single power-law LF jPaigne et al.l 
bOOd ). Figure O illustrates, in the case SFR3-I- Amati-like 
relation, the position of the best models in the parameter 
space of the LF. As can be seen, the low-luminosity slope is 
strongly constrained to be small, 5i <1, with a mean value 
S\ ~ 0.6, while the high-luminosity slope is larger, 82 > 1.4, 
with a mean value &2 ~ 1.7. The break luminosity (right 
panel) is not so well constrained with best models having 
Lb ==; 4 X 10^°-6 X 10^^ erg s"\ Figure g] compares the fit of 
the data points (log — log P diagram and distribution) 
with the best model obtained either with a power-law or a 
broken power-law LF. Both models are in good agreement 
with data, without a preference for one or the other. This is 
also indicated by the value of the reduced minimum in 
both cases : 1.4 (power-law) and 1.3 (broken power-law) for 
37 degrees of freedom. 

Figure [S] shows - for the best model - the LF as well as 
the luminosity distribution of bursts detected by BATSE, 
IIETE2 and SWIFT. It appears that the high-luminosity 
branch above Lb is extremely close to the best model single 
power- law LF (thin line). Below a few 10*^ erg s~^ (cor- 
responding to the lowest values of the constrast, k < 1.2), 
the fraction of detected GRBs is extremely low (less than 
10~* of the total). The two models (power-law vs broken 
power- law) differ in the lO^'^-lO'^^ erg s ^ range, where 



the fraction of detected bursts is still small (less than 20% 
of the total) with therefore little effect on the observable 
quantities. 

These results indicate that present data are not 
sufficient to distinguish between a power-law and a broken 
power-law LF. Both models can provide equally good fits 
to the observations. It is however interesting that a broken 
power-law remains allowed, as there are good theoretical 
arguments in favor of such a shape (see Sect. 2). Table [T] 
show that the properties of the broken power-law LF 
remain very stable as long as Lmin is kept to a low value 
(Lmin 10** erg s~^): the position of the break and espe- 
cially the values of the two slopes are not changing much, 
even for different GRB rates (SFRi and SFR2 have also 
been tested). It seems however that the broken power-law 
LF is more sensitive to the assumptions on the spectral 
parameters. In the case where the spectral properties are 
not correlated to the luminosity (log-normal peak energy 
distribution), the low- luminosity slope is not too different 
from the "Amati-like relation" case {Si ~ 0.7 instead of 
0.6), but the break luminosity is larger (Lb — 10^^"^'^ 
instead of io^°~^^ erg s~^) and the high-luminosity branch 
is steeper (82 — 1.8 — 2.4 instead of 1.7). 

These results can be partially compared to other 
studies. Based on an analys is of the BATSE log TV — log P 
diagram, IStern et al.l (|2002l ) have tested several shapes of 
GRB LFs, including a power-law and a broken power-law. 
Their assumptions concerning GRB spectra are different 
from those chosen in the present paper but are very close 
to our "log-normal peak energy distribution" scenario. For 
a GRB rate similar to our SFR3, they find 5i ~ 1.3 — 1.6, 
^2 > 3 and Lb ~ lO'^^ - 6 x 10^^ erg s"\ Whereas we are in 
reasonable agreement for the high luminosity slope, there 
is a large discrepancy for the low-luminosity branch, which 
is much steeper in their study. To understand the origin of 
this discrepancy, we made an additional simulation where 
we force 5i = 1.3 and let S2 and Lb free. We find that a 
good fit to the log A'^ — log P diagram can be found but 
that the peak energy distribution is not reproduced, which 
stresse s the importance of this additional constrain in our 
study. iFirmani et al.l l|2004l ) have presented a set of Monte 
Carlo simulations with assumptions on the intrinsic GRB 
properties that are very similar to ours but different obser- 
vational constraints, as they fit the distribution of BATSE 
pseudo-reds hifts obtained from the lu minosity-variability 
correlation (|Fenimore fc Ramirez-Rui3 2000). For the case 
with an intrinsic correlation between the luminosity and 
the spectral properties (Amati-like case) they find a break 
at Lb — 3 X 10^^ erg s~^ with the low and high-luminosity 
slopes equal to (Si ~ 0.8 and ^2 — 2.1. Taking into account 
the differences in the two approaches, the agreement 
between this study and our result is satisfactory , espe - 
cially for the slopes. More recently, iGuetta et al.l l|2005h : 
iGuetta fc PiranI l|2007l ) have also studied the GRB rate and 
LF using the recent results from SWIFT. Their analysis 
is based on the use of the log N — log P diagram only and 
they assume a very simplified GRB spectral shape, that is a 
power-law spectrum of photon index —1.6. For SFR3, they 
find a break luminosity Lb — 4 x 10^^ erg.s"^ and low and 
high luminosity slopes Si ~ 0.1 and S2 — 2, with large error 
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Figure 3. Parameter space (SFR3, Amati-like relation) : we plot the location of the best models (1 tr level) for a broken power-law 
LF with fixed minimum and maximum luminosities Lmin = 10'*^ erg s~^ and Lmax = lO^'^'^ erg (big dots). The range of parameters 
explored in the Monte Carlo is indicated with a box. In the right panel, the best models for a single power- law LF llDaigne et al.ll2006l) 
are also plotted with small dots (in this case the x-axis stands for the minimum luminosity and the y-axis for the slope). 




Figure 4. Best model (SFR3, Amati-like relation). Left: the simulated logAf — logP diagram of BATSE is plotted as well 
as BATSE data l lStern et al.l 12002 ); rig ht: the simulated peak energy distribution of bright BATSE bursts is plotted as well as the 
observed distribution dPreeceet"^ ] |2000h . In both panels, the best model for a broken power-law LF with fixed minimum and maximum 

and Lmax = 10®^'^ erg s~^ is plotted in thick line. For comparison the best model for a single power-law 
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LF obtained in lDaigne et al.l 1I2OO6I ) is plotted in thin line 
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Table 1. Best models; parameters. 



SFR 


log -Lmin 


log Ly~, 


log Lmax 


<5i 


^ 1 




log k 


log pQ 




(erg s-i) 


(erg s-^) 


(erg s-l) 






(keV) 




(GRB Gpc-3 yr-i) 








Amati-like relation E„ 


OC LO.43 








3^ 


50.3 ±0.7 




53.5 ± 0.4 


1.54 ±0.18 






-6.0 ±0.2 


-0.8 ±0.2 


r 


45 


50.4 ±0.5 


53.5 


0.65 ±0.22 


1.71 ±0.07 




-5.2 ±0.3 


0.0 ±0.3 


2 


45 


50.5 ± 0.4 


53.5 


0.62 ±0.20 


1.71 ±0.09 




—5.3 ± 0.2 


—0.1 ± 0.2 


3 


45 


51.2 ±0.6 


53.5 


0.60 ±0.22 


1.71 ±0.22 




-5.9 ±0.2 


-0.7 ±0.2 


3 


46 


51.1 ±0.7 


53.5 


0.64 ±0.20 


1.70 ±0.29 




-5.9 ±0.2 


-0.7 ±0.2 


3 


47 


51.2 ±0.6 


53.5 


0.69 ±0.24 


1.70 ±0.30 




-5.8 ±0.2 


-0.6 ±0.2 


3 


48 


51.3 ±0.5 


53.5 


0.74 ±0.27 


1.75 ±0.38 




-5.8 ±0.2 


-0.6 ±0.2 


3 


49 


51.6±1.1 


53.5 


1.02 ±0.45 


1.98 ±0.72 




-5.7 ±0.2 


-0.5 ±0.2 


3 


50 


52.0 ±1.2 


53.5 


1.47 ±0.65 


2.05 ±0.86 




-5.9 ±0.1 


-0.7±0.1 






log- 


-normal peak energy distribution 








3°^ 


50.5 ± 1.3 




53.7 ±0.9 


1.52 ±0.48 




2.79 ± 0.08 


-6.2 ±0.2 


-1.0 ±0.2 


1 


45 


51.2 ±0.4 


53.5 


0.67 ±0.23 


1.80 ±0.16 


2.80 ±0.11 


-5.5 ±0.2 


-0.3 ±0.2 


2 


45 


51.5 ±0.5 


53.5 


0.71 ±0.21 


2.04 ±0.57 


2.75 ±0.09 


-5.6 ±0.3 


-0.5 ±0.3 


3 


45 


52.1 ±0.5 


53.5 


0.66 ±0.20 


2.37 ±0.71 


2.80 ±0.08 


-6.1 ±0.2 


-0.9 ±0.2 



Preferred model in the case of a single power- law LF toaiene et al.ll2006t) . 



bars. There is therefore a discrepancy for the low-luminosity 
slope which is flatter in their case. We beheve however that 
the use of a more realistic spectral shape together with a 
constraint on the peak energy distribution leads in our case 
to a better determined LF at low luminosity. 

The two important results from this new set of Monte 
Carlo simulations are that (1) a broken power-law is 
compatible with present data but is not preferred compared 
to a single power-law (equally good fits of the observa- 
tions); (2) if the LF is indeed a broken power-law, the 
low-luminosity slope is constrained to be Si ~ 0.6 ± 0.2 
(Amati-like relation) or Si ~ 0.7 ± 0.2 (log-normal peak 
energy distribution), i.e. compatible with the prediction of 
the internal shock model derived in Sect. 2. 



3.3 The rate of underluminous GRBs 

A by-product of this study is an estimate of the local GRB 
rate, which is given in table [T] and is typically, in the SFR3 
scenario, 0.1 - 0.3 GRB Gpc"^ yr"^ in the Amati-like case, 
and 0.08 — 0.2 GRB Gpc~^ yr~^ for a log-normal peak en- 
ergy dis tribution. This i s in g ood agreement with the re- 
sults of iGuetta fc PiranI l|2007f ) . Despite the fact that our 
broken power-law LF can in principle extend to very low- 
luminosity, such a low local rate corresponds to less that 
10"'^ GRB yr~^ within 100 Mpc, which cannot explain the 
observation of GRB 9 80425 at z^O.008 or GRB 060 618 at 
z=0.03. As shown in iDaigne fc MochkovitchI l|2007l ). such 
underluminous bursts are well explained in the framework of 
the internal shock model by mildly relativistic / mildly ener- 
getic outflows. This would then indicate that the collapsing 
stars capable to generate such outflows, less extreme than 
those required to produce standard GRBs, are very numer- 
ous and should then produce an additional component in 
the LF, dominant at very low luminosity. This new branch 
cannot be the simple continuation of the LF of standard 



GRBs derived in this paper. Recently iLiang et all (|2007l ) 
have studied such a two component LF and found that the 
local rate corresponding to the low-luminosity component 
has to be several orders of magnitude above that of stan- 
dard GRBs but can still represent only a fraction of all type 
Ib/c supernovae. 



4 CONCLUSION 

We have demonstrated that in the framework of the internal 
shock model, a two branch LF is naturally expected, with 
a predicted low-luminosity branch which is a power-law of 
slope close to —0.5. This result is robust as long as the 
central engine responsible for GRBs is capable to produce 
a broad diversity of outflows, from highly variable to very 
smooth. 

Using a set of Monte Carlo simulations, we have then 
shown that current observations (log A'^ — log P diagram, 
peak energy distribution, fraction of XRRs and XRFs) are 
compatible with a broken power-law LF but still do not 
exclude a single power-law distribution. The low-luminosity 
slope of the broken power-law is strongly constrained to 
be Si ~ 0.4 — 0.9, compatible with the prediction of the 
internal shock model. 

These results are encouraging but only preliminary. 
A better determination of the GRB LF would provide an 
interesting test of the internal shock model when the low 
luminosity branch becomes more easily accessible. This will 
however require the difficult task of detecting many bursts 
at the threshold of current instruments and measuring their 
redshift and spectral properties. 
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Figure 5. Luminosity function in the scenario SFR3 + 
Amati-like relation. The apparent LF, as well as the lumi- 
nosity distribution of bursts detected by BATSE, HETE2 and 
SWIFT are plotted in thick line for the best model using a bro- 
ken power-law with fixed minimum and maximum luminosities 
^min = 10*^ erg s^^ and Lmax = 10^^-^ erg . All other pa- 
rameters can be found in table [T] For comparison the bes t model 
for a single power-law LF obtained in lDaigne et alJ ll2006h is plot- 
ted in thin line. Despite the fact that this figure was ob- 
tained with a special run simulating 10'' GRBs with the 
best model parameters, the curves are very noisy at low 
luminosity, as these events are very rare. 
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